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ABSTRACT 

We consider the consequences of gravitational wave recoil for unified models of active galactic 
nuclei (AGNs) . Spatial oscillations of supermassive black holes (SMBHs) around the cores of galaxies 
following gravitational wave (GW) recoil imply that the SMBHs spend a significant fraction of time 
off-nucleus, at scales beyond that of the molecular obscuring torus. Assuming reasonable distributions 
of recoil velocities, we compute the off-core timescale of (intrinsically type-2) quasars. We find that 
roughly one-half of major mergers result in a SMBH being displaced beyond the torus for a time of 
10 7 ' 5 yr or more, comparable to quasar activity timescales. Since major mergers are most strongly 
affected by GW recoil, our results imply a deficiency of type 2 quasars in comparison to Seyfert 2 
galaxies. Other consequences of the recoil oscillations for the observable properties of AGNs are also 
discussed. 

Subject headings: galaxies: active - galaxies: evolution - quasars: general 



1. INTRODUCTION 

Gravitational waves emitted anisotropically during 
gravitational collapse carry away linear momentum. As 
a result the center of mass of the collapsing object re- 
coils (Peres 1962; Bekenstein 1973). Configurations of 
coalescing spinning black holes can result in recoil ve- 
locities of hundreds to thousands of kms -1 (e.g., Cam- 
panelli et al. 2007; Gonzalez et al. 2007; Herrmann et 
al. 2007; Pollney et al. 2007; Tichy & Marronetti 2007; 
Brugmann et al. 2008; Dain et al. 2008), scaling to a 
maximum of ~ 4000 kins -1 (Campanelli et al. 2007; 
Baker et al. 2008) for maximally spinning equal-mass 
binaries with anti-aligned spins in the orbital plane, and 
as large as ~ 10 4 kms -1 in hyperbolic encounters (Healy 
et al. 2008). Kicks large enough to remove SMBHs from 
galaxies have potentially far-reaching consequences for 
SMBH and galaxy assembly, and predict interstellar and 
intergalactic quasars (e.g., Madau et al. 2004; Merritt et 
al. 2004; Madau & Quataert 2004; Haiman 2004; Boylan- 
Kolchin et al. 2004; Libeskind et al. 2006; Loeb 2007; 
Schnittman 2007; Volonteri 2007; Gualandris & Merritt 
2008; Volonteri & Madau 2008, Kornreich & Lovelace 
2008). 

Komossa et al. (2008) reported the detection of a re- 
coil candidate with a projected kick velocity of 2650 km 
s" 1 . The quasar SDSSJ092712. 65+294344.0 shows three 
unusual emission-line systems, including a kincmatically 
offset broad-line region (BLR) and a system of atypically 
narrow emission lines which lack the usual ionization 
stratification - two key signatures of kicks. Apart from 
spectroscopic signatures (Merritt et al. 2006; Bonning et 
al. 2007), offset quasars could be detected for instance 
by their temporarily flaring accretion disks (Shields & 
Bonning 2008; Lippai et al. 2008; Schnittman & Kro- 
lik 2008), by tidal disruption flares from the bound (and 
unbound) population of stars and episodic fuelling from 



stellar mass loss (Komossa & Merritt 2008), and via the 
bound compact star cluster itself (Merritt et al. 2008, 
O'Leary & Loeb 2008). One key consequence of gravita- 
tional wave recoil is long-lasting oscillations of the SMBH 
about the galaxy core, implying that SMBHs may spend 
as long as 10 6-9 yrs off-nucleus with an amplitude of par- 
sees or kiloparsecs, depending on kick velocity and galaxy 
structural parameters (Merritt et al. 2004; Madau & 
Quataert 2004; Gualandris & Merritt 2008; MQ04 and 
GM08 hereafter). In this Letter, we consider conse- 
quences of these "recoil oscillations" for unified models 
of active galactic nuclei (AGNs). 

According to unified models, AGNs are intrinsically 
similar, but their appearance depends strongly on the 
line of sight of the observer toward the "central engine" 
(reviews by Antonucci 1993; Elitzur 2007). Along certain 
sightlines, a dusty torus consisting of molecular clouds 
blocks the observer's view, hiding some core components, 
especially the BLR, which therefore is only directly visi- 
ble in "type 1" AGNs. The unified model has been very 
successful in explaining observed properties of AGNs in- 
cluding the presence of hidden BLRs detected in polar- 
ized light (e.g., Antonucci & Miller 1985; Zakamska et 
al. 2005), and has been corroborated by recent imaging 
and spectroscopy of the torus (e.g., Jaffe et al. 2004; 
Siebenmorgen et al. 2005). 

If SMBHs and the BLRs bound to them spend a signif- 
icant time displaced from the nucleus, this will have pro- 
found consequences for obscuration-based unified models 
of AGNs. This statement holds true whether the ob- 
scuration originates in parsec-scale molecular tori (An- 
tonucci 1993), or compact star-forming regions (Leven- 
son et al. 2001), or is due to absorption associated with 
the host galaxy itself, which may work on typical scales 
of ^100 pc (Maiolino & Rieke 1995). Here we argue 
that a significant fraction of the quasar population is ex- 
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pected to be in a regime such that the SMBH and the 
BLR bound to it is displaced beyond the obscuring re- 
gion, implying a deficiency of type 2 (obscured) AGNs 
among quasars. 

2. RELEVANT SPATIAL SCALES IN THE AGN CORE 

We first compare three characteristic spatial scales: the 
kick radius, the BLR size, and the torus size. After the 
kick, matter remains bound to the SMBH within a region 
whose radius is given by 



GM bh nAO f M BH 
rk = 5 — w 0.43 



VlO 8 M y Vl0 3 km s 



pc 



(1) 

where is the kick velocity (e.g., Merritt et al. 2006). 
The size tblr of the BLR of AGNs has been determined 
from reverberation mapping (Peterson 2007), and scales 
with AGN luminosity as 
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(Kaspi et al. 2005) where the luminosity at 5100A, 
AL A (5100A) ~ 0.1 L and L is the AGN bolometric lu- 
minosity. The size of the molecular torus is still poorly 
constrained from observations. Recent measurements of 
dusty gas in the Seyfert galaxies NGC 1068 and Circinus 
suggest an extent of 2-3 pc (Jaffe et al. 2004; Davies et 
al. 2007; Tristram et al. 2007). It is reasonable to as- 
sume that the inner edge of the torus, r tor ,in, is beyond 
the dust sublimation radius (e.g., Netzer 1990; Nenkova 
et al. 2008) which is given by 

L \ 1/2 /1500K\ 2 6 
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with an outer radius likely not much larger than ~ 
20?"tor,in (Elitzur 2007), and where T su b is the dust subli- 
mation temperature. Comparison of these three relations 
shows that a large fraction of the BLR remains bound 
to the recoiling hole, while structures of the size of the 
torus or larger will typically be left behind. Oscillation 
amplitudes of <; 10-20 pc will therefore move the SMBH 
beyond the torus scale, except for the highest SMBH 
masses, where part of the torus will remain bound. 

3. OSCILLATION AMPLITUDES AND OSCILLATION 
DURATIONS IN QUASARS 

We base our discussion on the A-body simulations of 
GM08, who computed SMBH trajectories after GW re- 
coil. Key model parameters are the SMBH mass Mbh 
and galaxy mass M gal = 10 3 M B h> the kick velocity v^, 
and the galaxy structural parameters. We concentrate 
here on the mass range of SMBHs that is typical for the 
bulk of quasars (M B h « 10 8 - few 10 9 M o ). Massive in- 
active elliptical galaxies and host galaxies of luminous 
quasars typically have effective radii R e on the order of 
a few to ~ 10 kpc. Here, we adopt the relation between 
R e and M ga i defined by the nine luminous (— 22 < My < 
—24) quasar host galaxies in the sample of Wolf & Shcinis 
(2008), i.e. log(i? c /kpc) * -5.61 + 0.55 log(M gal /M Q ). 
The models of GM08 on which we base our discussion had 
post-kick core radii r c consistent with the range defined 
by the brightest E-galaxies with resolved cores (Ferrarese 
et al. 2006), i.e. 0.01 < r c /R c < 0.03. We focus here on 



the spherical "Al" galaxy models from GM08; scaling of 
the A-body results to physical units was done following 
their equation (4). Dark matter halos were ignored. 

Figure la shows trajectories of kicked SMBHs scaled 
to a galaxy with SMBH mass 5 x 10 8 M Q . Kick ve- 
locities were Vk/v csc = 0.3,0.5,0.7,0.9, corresponding to 
Vk = (360,590,830,1070) km s _1 for a central escape 
velocity v csc of 1185 km s _1 of our galaxy model. As 
discussed by GM08, SMBH oscillations continue well be- 
yond the time ("Phase I") that would be predicted by 
applying Chandrasckhar's dynamical friction formula as- 
suming a fixed galaxy core (e.g. MQ04; Blecha & Loeb 
2008). The time of onset of these long-term, or "Phase 
II" , oscillations is indicated by the open circles in Fig- 
ure la. Including the effect of the Phase II oscillations, 
Figure la shows that the SMBH's motion persists for 
more than ~ 10 7 yr if t>k > 450 km s _1 . 

4. COMPARISON WITH OBSERVATIONS 

4.1. Unified models and typel/type2 fractions 

About 70% of the nearby Seyfert galaxies are type 2 
(e.g., Schmitt et al. 2001), i.e. they lack a BLR in their 
(unpolarized) optical spectra. The fraction of the high- 
luminosity equivalents, type 2 quasars, is less well known 
and still subject to a number of selection effects (e.g. 
Halpern & Moran 1998; Reyes et al. 2008), even though 
studies generally indicate a deficiency of type 2 quasars. 
While in past studies type 2 quasars were observation- 
ally very rare or absent, significant numbers have been 
found in recent X-ray and optical surveys (e.g., Norman 
et al. 2002; Zakamska et al. 2005; Brusa et al. 2007). 
There is a systematic trend such that the fraction of type 
2 sources, or equivalcntly the amount of X-ray absorp- 
tion, decreases with increasing source luminosity (e.g., 
Simpson 2005; Barger et al. 2005; Reyes et al. 2008; 
Hasinger 2008), while the simplest possible version of 
the unified model would imply a constant typel/type2 
ratio. Models have been proposed which account for 
this luminosity dependence, e.g. by invoking changes 
in the properties of the obscurer as L increases (e.g., 
Nenkova et al. 2008; Ballantyne 2008). Recoil oscilla- 
tions inevitably affect the numbers of obscured versus 
unobscured sources and therefore have potentially pro- 
found implications for unified models. How do they af- 
fect the ratio of typel/type2 quasars in comparison to 
the number of typel/type2 Seyfert galaxies ? There is 
increasing evidence that quasar activity is powered by 
major mergers while Seyfert activity may have other trig- 
gers including bars, minor mergers, or random accretion 
of molecular clouds (e.g., Sanders et al. 1988; Urrutia 
et al. 2008; Hopkins et al. 2008; Hasinger 2008, and 
references therein). Recoil oscillations have highest am- 
plitudes in major mergers, and we may therefore expect 
that the frequency and properties of (type 2) quasars are 
strongly affected by recoil oscillations. Can this explain 
the relative scarcity of type 2 quasars in comparison with 
type 2 Seyferts, and the more general trend that type2- 
ness (X-ray absorption) decreases with luminosity? 

4.2. Rate estimates 

Since a large fraction of quasars are believed to be trig- 
gered by major mergers, in order to estimate the fraction 
of quasars that occur at or above a given kick velocity, 
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we carried out rate estimates relevant for major mergers 
with random spin distributions (Campanelli et al. 2007; 
Schnittman & Buanano 2007; Baker et al. 2008; note 
that the actual kick velocities could be smaller in gas- 
rich systems if the mechanism discussed by Bogdanovic 
et al. 2007 is at work, but see Sect. 3 of Schnittman & 
Krolik 2008). We first updated these previous estimates, 
based on the recent kick formula of Baker et al. (2008; es- 
sentially identical results were obtained using the Lousto 
& Zlochower 2008 version of the kick formula), and as- 
suming random orientations of the spin vectors of both 
SMBHs and a distribution of SMBH mass ratios in the 
range 0.3 < q < 1.0 where q = 7712/7711,7712 < mi, rele- 
vant for major mergers. SMBH spins were drawn from 
a distribution such that a\ < a < 02 with a = S/m 2 the 
dimensionless spin, with ai — 0.5, a 2 = 0.9. The distri- 
butions of m and a were assumed to be uniform in the 
logarithm between these limits. For this base model, we 
find that about 50% of major mergers have kick veloci- 
ties above 500 km/s. At or above these kick velocities, 
total oscillation timescales start to be on the order of 
quasar lifetimes, which are about 10 7-8 yr (e.g., Kauff- 
mann & Haehnelt 2000; Yu & Tremaine 2002; Hopkins 
et al. 2005; see review by Martini 2004). 

How much of the total quasar population is ultimately 
affected by recoil oscillations then depends on (1) the 
time t v is the kicked SMBH+BLR spends beyond the ob- 
scuring torus (i.e., at distances greater than r tor us ~ 
^tor,out from the galaxy center (eqn. [3])) and therefore 
appears as "type 1" rather than "type 2" , in comparison 
to (2) the total quasar lifetime. We computed t V [ S on a 
grid in r tor us for each of the nine A-body trajectories in 
GM08 (0.1 < Vk/fcsc < 0.9). Given arbitrary values of 
^k, -Mg a i and r tor us, the value of t v - ls was then computed 
via interpolation between the nine discrete kick velocities 
of GM08. 

Figure lb shows, for our base model, the fraction 
/ of kicked SMBHs that spend more than 3 x 10 7 yr 
at r > r m ; n , for a range of (Af ga i,r m i n ) values. For 
Tmhi ~ ftor,out, this fraction is ~ 0.5 ± 0.1 with a weak 
dependence on galaxy mass. We conclude that a signifi- 
cant fraction of mergers would result in recoiling SMBHs 
that spend of order a quasar lifetime above the obscuring 
torus and would therefore appear as type 1 quasars even 
if "intrinsically" type 2. We examined the robustness 
of these results under different assumptions about the 
pre-recoil distributions of spins and masses. Maximally 
spinning, equal-mass SMBH mergers affect a fraction of 
up to / ~ 0.75 of the population (same assumptions as 
above), while major mergers with intermediate SMBH 
spins of a = 0.3 imply / ~ 0.05 — 0.15. 1 

If we assume that the typel /type2 fraction of quasars 
is intrinsically (i.e., in the absence of recoils) the same 
as in Seyferts (~70% type 2, ^30% type 1), we predict 
that a fraction ~ 50% of these type2s will appear in- 
stead as type Is at any given time. Several factors affect 
these estimates: idealizations in the galaxy models used 
in the -ZV-body simulations, uncertainties in the distri- 
bution of masses and spins as discussed above, uncer- 

1 We note that galaxy triaxiality, not included in the current 
models, might extend the oscillation timescales by a factor of sev- 
eral (e.g. Vicari et al. 2007; see also MQ04 and Merritt et al. 2004); 
inhomogeneities of the host galaxy after merging would likely have 
the same effect. 



tainties especially in the thickness of the torus in the 
most massive galaxies, and various selection effects in 
the measurements of typel/type2 ratios in dependence 
of SMBH mass. We also recall that if equations (1) and 
(3) strictly hold, in the most massive quasars a fraction 
of the torus will remain bound to the recoiling SMBH so 
absorption/extinction would not fall to zero. Finally, we 
note that if the quasars' total lifetime is actually com- 
posed of several shorter merger episodes on the order 
of 10 6 yr each, even recoil oscillations with velocities as 
small as ~ 200 km s" 1 would affect a large fraction of 
the quasar population. 

4.3. Implications 

So far, we have distinguished between Seyfert galax- 
ies (low-mass SMBHs) and quasars (high-mass SMBHs) . 
Can we also reproduce the observed trend (Sect. 4.1) 
that obscuration fraction decreases with quasar luminos- 
ity (i.e., mass)? In our base model, the dependence of 
mean oscillation timescale on galaxy mass is weak. How- 
ever, since the likelihood of a major merger (as opposed 
to other types of fuelling) is believed to increase strongly 
with AGN luminosity (galaxy mass), the observed trend 
should arise naturally, since the most luminous AGNs 
are increasingly likely to be triggered by major mergers. 

If a fraction of all quasars is recoiling at any given time, 
we should see the corresponding BLR emission-line veloc- 
ity shifts w bs in a fraction of the type 1 quasars. Bonning 
et al. (2007; B07 hereafter) set limits on the fraction of 
emission line velocity shifts observed in a sample of SDSS 
quasars. In our model, the majority (> 75%) of kicked 
SMBHs remain bound to the galaxy and so their veloc- 
ities quickly drop below the initial value of v^; most of 
the time they would therefore be observed with a velocity 
that is much smaller than v^. 

We carried out Monte-Carlo simulations to check the 
consistency of our recoil model with the B07 limits. Our 
simulations were designed to crudely mimic the prop- 
erties of a sample of quasars selected to exhibit both 
broad and narrow emission lines, as in the B07 sample. 
We adopted a uniform logarithmic distribution of galaxy 
masses, 11 < log 10 (Mg a i/M Q ) < 12.5, and a fixed quasar 
lifetime of t qso — 3 x 10 7 yr; kicks were assumed to have 
occurred at times that were distributed uniformly and 
randomly between the epoch of observation and a time 
i qso earlier. The position and velocity of the SMBH at 
the time of observation was extracted from the appropri- 
ate iV-body model after scaling to physical units using 
the galaxy mass. If the distance of the SMBH from the 
galaxy center exceeded the outer torus radius, its radial 
velocity u b s was added to the cumulative distribution 
(assuming a random direction for the recoil); this veloc- 
ity was identified with the measured velocity offset of the 
BLR gas in the B07 galaxies. A fraction 30% of SMBHs 
(the "true" type 1 population) with r < r tor! out were as- 
sumed to have visible BLRs and so were included in the 
accounting. Figure lc shows that the predicted fraction 
of objects with large (<; 100 km s _1 ) velocity shifts is 
~ an order of magnitude smaller than would be inferred 
from the unmodified distribution of kick velocities; this is 
due to the deceleration that occurs as the SMBH moves 
through the galaxy. 

B07 found a maximal fraction of /5oo = 0.04 quasars 
with velocity shifts above v b s = 500 km s , and a frac- 
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tion /iooo=0-0035 with shifts above u bs = 1000 km s . 
These limits are consistent with our baseline model for 
the kick distribution, particularly if we impose an upper 
limit to the galaxy/SMBH separation at which a recoil- 
ing SMBH would be spectroscopically identified with its 
host galaxy (shown as the dotted lines in Figure lc, as- 
suming r max = 10 kpc). 

Finally, we note that recoil oscillations will also have 
a number of other observable consequences. They will 
affect the X-ray background and its modeling since a 
fraction of sources will be unobscured at any given time. 
In particular, small amplitude oscillations of the order 
the torus size will affect the ratio of Compton-thin to 
Compton-thick sources, and could lead to measurable 
variability in the absorption and extinction of AGN spec- 
tra once the recoiling SMBH passes the individual clouds 
making up the torus. A number of interesting effects are 
related to the torus itself: (1) A recoiling SMBH with 
a bound gas disk that passes through the dense torus 
(rather than moving perpendicular to it) might cause lo- 
cal shocks, heating, and temporary X-ray emission. (2) 
During the long-lived "Phase II" oscillations, when the 
SMBH oscillation amplitude is on the torus scale, the 
SMBH might efficiently accrete from the dense molec- 
ular gas at each turning point, causing repeated flares 
of radiation. Such flares would locally destroy the dust, 
while photoionization of the dense surrounding gas would 
produce a strong emission-line response which would not 

2 While strong emission-line variability in response to an X-ray 
flare has recently been observed, this event is more likely inter- 
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only help in identifying recoils but could also be used 
as a new probe of the properties of the torus itself. 2 
(3) Torus radii are roughly equal to SMBH gravitational 
influence radii, so ejection of the SMBH might lead to 
temporary expansion of the torus since the mass hold- 
ing it in place is suddenly removed. Isotropic expansion 
would not affect the column density along the line of 
sight, seen from the very center. These effects will be 
addressed in more detail in forthcoming work. 

In summary, we have shown that timescales of recoil 
oscillations are in an interesting regime where they can 
potentially affect a significant fraction of the quasar pop- 
ulation. Details of the predictions still depend on uncer- 
tainties in the observed numbers of AGN types, quasar 
lifetimes, torus properties, structural parameters of lumi- 
nous quasar host galaxies, and SMBH spin distributions 
on the one hand, and on modelling recoil in realistic non- 
spherical or non-axially-symmetric galaxies on the other 
hand. Knowledge of recoil oscillation timescales and am- 
plitudes is also critical for modelling AGN evolution, de- 
lays between starburst and AGN activity after merger, 
and the cosmic X-ray background. 
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time [yr] log M BH [M sun ] V [km/s] 

Fig. 1. — (a) Trajectories of kicked SMBHs in N-body models of a galaxy (based on GM08) with a SMBH mass of M B h = 5 x 10 8 M q , 
galaxy mass M ga \ = 5 X 10 Mq, and effective radius of 7 kpc. Kick velocities were (360, 590, 830, 1070) km s" 1 (blue, red, green, black). 
Open circles indicate the approximate time at which the SMBH would come to rest after "Phase I" (see text), (b) Fraction of kicked 
SMBHs, in a galaxy with SMBH mass Mbh, that remain above the torus, of radius r m ; n , for a time of 3 X 10 7 yr or longer. Thick (green) 
contours are based on the full ./V-body trajectories; thin (white) contours include only "Phase I." Gray-scale density is proportional to 
the fraction of kicks that result in escape, from ~ 0.3 on the left to ~ 0.05 on the right. Kicks were generated assuming "log-uniform" 
distributions of SMBH masses and spins (see text) with 0.3 < m2/mi < 1 and 0.5 < < 0.9. The red line shows the approximate 

outer radius expected for the obscuring torus (based on eqn (3), and converting Mbh to L assuming accretion at 0.1L e( jd)- (c) Velocity 
distributions that would be observed in a representative sample of quasars, assuming that kicks occur randomly in time; additional features 
of the model are described in the text. The black line uses the baseline model of SMBH masses and spins in the pre-recoil binary (the same 
model used in Fig. lb). Red line: like the base model, but a = 0.3 for both SMBHs. Blue line: equal-mass SMBHs with maximal spins. 
The dotted lines exclude SMBHs that are more than 10 kpc from the galaxy center at the moment of observation. The gray dashed line 
shows the input kick distribution for the baseline model, unmodified by motion through the galaxy. 



